Crossovers in spin-boson and central spin models 
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We discuss how the crossovers in models like spin-boson model are changed by adding the coupling 
of the central spin to localised modes- the latter modelled as a 'spin bath'. These modes contain 
most of the environmental entropy and energy at low T in solid-state systems. We find that the low 
T crossover between oscillator bath and spin bath dominated decoherence, occurring as one reduces 
the energy scale of the central spin, is characterised by very low decoherence- we show how this 
works out in practise in magnetic insulators. We then reconsider the standard quantum-classical 
crossover in the dynamics of a tunneling system, including both spin and oscillator baths. It is 
found that the general effect of the spin bath is to broaden the crossover in temperature between 
the quantum and classical activated regimes. The example of tunneling nanomagnets is used to 
illustrate this. 
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Is INTRODUCTION 

In the book of Weiss on quantum dissipative phenom- 
ena one finds a very nice summary of results on the 
crossover between quantum tunneling and classical ac- 
tivation for a single tunneling coordinate coupled to a 
bath of oscillators (see Chapters 10-17, particularly 14 
and 16) . This kind of problem has a long and interest- 
ing history, beginning with work of Kramers |2( in 1940. 
The oscillator bath models assume that each bath mode 
is weakly perturbed, and then the description of the bath 
by oscillators is well known to be correct. Many physi- 
cal j^ystems are very accurately described by such models 
ULId)Q>la> an d they are central to much of reaction rate 
chemistry as well. Typically one studies either a particle 
tunneling from a trapped state to an open continuum of 
states (the dissipative tunneling problem), or a double- 
well system in which a particle has to go from one well 
to another (the dissipative 2- well problem). One has a 
range of temperatures in which both activation and tun- 
neling processes are important. Both the width of the 
crossover regime and the detailed dependence of transi- 
tion rates, as a function of temperature and applied bias, 
are of interest 0, • In the 2 well problem, the 'quan- 
tum limit', where only the 2 lowest levels of the 2- well 
system are relevant (assuming a weak bias between the 
wells), has been studied very extensively. This is the 
'spin-boson model', in which a 2-level system couples to 
the oscillator environment. 

Another interesting application of the spin-boson 
model is to the problem of qubits in quantum information 
processing (QUIP). The central issue here is the study of 
decoherence in the dynamics of the qubit, and how it 
depends on both simple things like applied fields, tem- 
perature, etc., and in a more complex way on the de- 



tailed nature of the bath, and its coupling to the qubit. 
It turns out that at the low temperatures that are ap- 
propriate for QUIP, or for any other large scale quantum 
coherence, the oscillator bath models are no longer ade- 
quate to describe all the physics. In many systems the 
decoherence is controlled largely by the coupling to lo- 
calised modes, such as defects, tunneling charges, para- 
magnetic spins, or nuclear spins, and this environment 
of localised modes cannot in general be modelled by 
oscillators- it can however be described as a set of spins 
(the 'spin bath'). There is now extensive experimental 
evidence for the key role of such modes in experiments on 
Cooper pair box qubits Jl , SQUID qubits 0, 0] , and in 
molecular magnets |llL Il2j , but the importance of these 
modes is already rather obvious just from an estimation 
of their coupling to these systems. There have been a fair 
number of theoretical studies of spin bath environments. 
Early partial studies, in various contexts, include refs. 
[l3L H4I Hfl Hif : later work has concentrated on appli- 
cation to coherence and relaxation in tunneling systems 
(see, eg., refs. E3 E3), to decoherence in 

mesoscopic conductors I22L 123) . superconducting qubits 
0,I3> anc I nanomagnets |25|7lt is clear that in the limit 
of weak coupling to the spin bath, it should be possible 
to map to an oscillator bath- studies in this limit appear 
in, eg., refs. 0, LU4 |2fJ, |23 

One can think of this breakdown of the oscillator bath 
model in several ways. One is dealing here with a break- 
down of the assumption of weak system-bath couplings, 
and a corresponding breakdown of linear response in the 
behaviour of the bath dynamics (for more detailed dis- 
cussion of this see ref. |jfl, and refs. therein, and also 
the nice short summary by Weiss, in ref. 0, pp. 49- 
52). We emphasize that the coupling to localised bath 
modes is almost always weak compared to the tunneling 
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barrier energy Eb , or to the energy lo corresponding to 
small oscillation energies in the potential wells- so these 
modes are usually invisible in ordinary tunneling exper- 
iments. However, the energy scale of the localised bath 
modes, and their coupling to the central qubit coordi- 
nate, is often not small compared to the exponentially 
smaller tunneling energy A Q . In the qubit regime it is 
the comparison with A Q that counts, particularly for de- 
coherence (for a more precise discussion see below). 

We are thus left with an interesting problem. What 
is the combined effect of spin bath and oscillator bath 
modes on the dynamics of the central system? In partic- 
ular, how is decoherence affected by these two, and how 
does the oscillator bath take over from the spin bath as 
one goes to higher temperatures, or increases A Q ? Es- 
sentially one has to reconsider the whole question of the 
crossover between quantum and classical regimes when 
both baths are included. 

In this paper we give a progress report on these ques- 
tions for 2 kinds of crossover, viz: 

(i) The crossover between spin bath controlled decoher- 
ence, which dominates when A Q is small, and the oscilla- 
tor bath-controlled decoherence, which dominates when 
A c is large (in both cases, assuming low temperature). 
The most interesting behaviour is in the crossover regime 
itself, when the decoherence goes to a minimum. Thus by 
raising A Q one can go from an incoherent regime, through 
a regime of coherent qubit dynamics, and then back to 
incoherent tunneling. To illustrate the idea we show how 
the general idea works for nanomagnets coupled to nu- 
clear spins and phonons- The detailed application to spe- 
cific magnetic and superconducting systems is discussed 
elsewhere [25tl28j . 

(ii) We look at how the spin bath influences the 
crossover between the quantum tunneling and classical 
thermally activated regime. This also involves a crossover 
between spin bath and oscillator bath environments. 
Given the complexity of this crossover, we do not attempt 
any complete discussion, but instead make some qual- 
itative remarks on the physics, and then present some 
results for magnetic insulators (again involving phonons 
and nuclear spins) . For related work one may go to a 
series of papers |2^| on the application to ensembles of 
tunneling magnetic molecules. 



2: CROSSOVERS FROM COHERENCE TO 
INCOHERENCE 

The spin-boson model has a control parameter A Q (the 
operating frequency of the qubit); and we consider here 
the crossover between the small A G regime, where deco- 
herence is controlled by the spin bath, and the large A Q 
regime, where it is controlled by the oscillator bath. The 
interesting thing is that in the crossover between these 2 
regimes lies a 'dead zone' where decoherence can be very 



low. This 'coherence window' will be very important for 
solid-state based quantum information processing. 



2a: QUBIT COUPLED TO OSCILLATOR AND 
SPIN BATHS 

We consider a 2-level system (a qubit) with the usual 
bare Hamiltonian 



H = A c 



(1) 



This is coupled to both spin and oscillator baths. The 
thermal energy fc^T, and the longitudinal and transverse 
field energies e ,A , are assumed to be much less than 
the energy gap E g to any higher levels of the system. In 
a magnetic qubit (eg., a magnetic molecule, or a rate 
earth ion), this 'spin gap' is typically 5 — 10 K, and 
in a superconducting qubit the corresponding Josephson 
plasma frequency depends strongly on the junction ge- 
ometry, and might be a little less. 

The baths themselves are assumed to have Hamiltoni- 
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in terms of a set of oscillators {x q } describing delocalised 
modes and a set of spins {o~k} describing localised modes 
(here for simplicity assumed to be a set of Pauli spin-1/2 
systems). We have written the set of 'fields' {hj,}, acting 
on the individual bath spins, in the form hj, = co^rhk, 
where rhk is a unit vector in the direction of the field. We 
assume that a UV cutoff Q a exists in these Hamiltonians, 
so that all spin and oscillator degrees of freedom have 
energy < Q a . The 2 baths are coupled to the central 
qubit via the following diagonal couplings: 



(4) 



where {4} are a set of unit vectors. There can also 
be non-diagonal couplings, ie., terms which operate only 
when the qubit is switching between the eigenstates | j) 
and | |) of f z . These are usually specified by modifying 
the form of the transverse term A t x in the bare qubit 
Hamiltonian. For the oscillator bath one adds a coupling 

m 



i 



-(r + c q x q + H.c.) 
and for the spin bath one makes the substitution 



Aofj. 
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A f + e' 



H.c. 



(5) 
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where the {n^} are unit vectors. 

We briefly note the important features of these inter- 
actions. First, we recall that the usual longitudinal spin- 

I —1/2 

boson couplings {cij} are typically ~ 0(N ), where 
N a is the number of oscillator degrees of freedom in the 
Hilbert space (defined by the UV cutoff Q ). On the 
other hand the spin bath couplings {wji} may have a 
quite different dependence- in magnetic qubit systems 
they are usually independent of N s , the number of spins, 
whereas in SQUID qubit systems one has uj\ ~ 0(N~ 2 ), 
at least in the simplest designs. For large N Q this means 
that the oscillator bath couplings are very weak- justify- 
ing the initial model of linear weak couplings. In the case 
of a SQUID qubit coupled to a spin bath one sees that it 
ought to be possible to map t he p roblem to a spin-boson 
model, and indeed one can |7l.l28|. However this is not an 
option for magnetic qubits- not only are the individual 
hyperfine couplings between the qubit and the nuclear 
spins independent of the number N s of nuclear spins, 
they are also large- in many cases w\ for a single nuclear 
spin can exceed A ! In this case we must deal directly 
with the spin bath, and give up any hope of mapping the 
problem to a spin-boson model. 

A second remark concerns the non-diagonal couplings. 
In cases where the diagonal couplings happen to be zero 
(which can happen under unusual circumstances) the 
non-diagonal couplings are the only remaining decoher- 
ence mechanism- this makes them very interesting for 
studies of decoherence (a point which has also been noted 
in recent discussions of superconducting qubits H3]). On 
the other hand when the diagonal couplings are non- 
zero, they usually dominate over the non-diagonal ones, 
at least when the qubit is modelling a tunneling solid- 
state system. It then follows that both qf /c|J and a% are 

small- in fact c^/cf ~ 0(A O /O O ) and ~ 0(w|/0 o ) 
(for more details on this see refs. pHli||). 

Finally, we note that the interactions between bath 
modes are treated differently in the oscillator and spin 
bath cases. In the oscillator bath case it is is usually 
argued that any weak anharmonic interactions have lit- 
tle relevance to the dissipation or decoherence caused 
by the bath- that information and energy are quickly 
transported away from the qubit, and so we can drop all 
reference to intra-bath interactions. In the case of the 
spin bath, however, it is clearly incorrect to drop such 
interactions- even though they are usually very small. 
This is because the spin bath describes local modes, 
which are not weakly coupled to the qubit- accordingly a 
large amount of energy and information can in principle 
be dumped into each mode and the V^cr^cr^, interaction 
is the only way this can be redistributed. Over long time 
scales non-linear effects become inevitable, and the size 
of becomes very important. We shall see how this 
works below. 



2b: DECOHERENCE RATES 

We define the decoherence dynamics for the qubit in a 
fairly standard way by assuming an initial state | t), 
and calculating the reduced density matrix as a function 
of time thereafter, once the spin and oscillator baths are 
integrated out. The general form of the result (assuming 
the bias e a — for simplicity) is 

2p(t) = 

/(l+ e - r i«cos(2A t)) ie~ v ^ sin(2A t) \ 

V -ie- T ^ sin(2A *) (1 - e- ri « cos(2A t)) / 

(7) 

Here 1^(4), with (i = 1,2, may have a complicated 
time dependence. If r M (i) — * = const., (so that the 
coherence decays exponentially in time) one can write 
l/r M = T M , following NMR terminology. In this case 
we say that, in this basis, the decay rate I/T2 of the 
off-diagonal matrix elements is the decoherence rate- it 
characterizes the rate at which interference between | f) 
and I I) states is lost. In other cases one can usually 
derive a characteristic timescale for the loss of coher- 
ence, and this is called the decoherence time. One also 
defines a dimensionless measure of the decoherence rate, 
given by 

7 = I/r A o (8) 

or its inverse, the 'decoherence quality factor', often de- 
fined as Q0 = n /jrf> = ttt^Aq. This Q-factor tells us 
roughly the number of coherent oscillations of the sys- 
tem before decoherence sets in. 

Here we first quickly recall the known results for deco- 
herence in this kind of problem. The dimensionless decay 
rate 7^ has the following contributions: 

(i) Oscillator bath contributions: The decoher- 
ence rates here depend on the Caldeira-Leggett spectral 
density . For the cases we are interested in one has the 
following results: 

(a) Phonon decoherence: this is relevant when we deal 
with spin-phonon coupling in magnetic insulators. Typ- 
ically one considers a spin S (representing a molecular 
spin or other nanomagnet), which truncates at low tem- 
perature to a magnetic qubit (when ksT <C fl , where fl 
is the spin gap to the higher electronic excitations). For 
decoherence the most important coupling to S is the non- 
diagonal coupling ~ ST2 |q| to acoustic phonons having 
a Debye energy 80 (for a simple derivation of this see ref. 
(lif. and for a thorough discussion of spin-phonon cou- 
plings see Standard spin-boson methods 0,0, 
applied to this coupling, give a contribution 7? to 7^ 
of perturbative (ie., golden rule) form; when the applied 
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bias is zero one gets (l8| : 

[{sn A f/e 4 D } coth(A /k B T) 



ph 



(9) 



which is very weak at low energies (ie., for A a <C #d). 
Although no qubit behaviour has yet been seen in mag- 
netic systems, there are extensive experimental results for 
the effect of phonons on the spin dynamics of magnetic 
molecules. 

(b) Electronic decoherence: This comes in when we 
need to analyse decoherence in SQUID flux qubit sjojjld 
El Hi] or Cooper pair box charge qubits [1 Hi EJ 
For example, in flux qubits tunneling between flux states 
±0 m , and with charging energy E c , one has a dimen- 
sionless coupling a — (Idxf^ujo/ E^Q^ 1 between SQUID 
flux and electronic bath, parametrised in terms of the 
SQUID Q-factor. In this case of Ohmic dissipation the 
decoherence rate is |l| 



7 ;(A D ) = -acoth(A /fc B T) 
where again we assume the system is in resonance. 



(10) 



(ii) Spin bath contributions: We write the contri- 
butions to the spin bath-induced decoherence in terms of 
the couplings introduced in @ and JBJ. It is useful to 
also introduce another quantity E Q which quantifies the 
total effect on a single qubit level of the coupling to the 
bath spins- we have 



E, 



(ii) 



so that E a is just the half- width of the Gaussian envelope 
of 2 Na spin bath states associated with each qubit state. 
This formula is easily generalised to include 'higher spin' 
bath spins (see below). In section 3 we say more about 
the structure of these bath spin states inside this 2 Ns -fold 
manifold. 

There are 3 spin bath contributions to the decoherence 
00: 

(a) Noise decoherence: spin diffusion inside the spin 
bath causes the longitudinal bias acting on f z to fluctu- 
ate over a range Tm in energy bias space, causing phase 
noise. There will be a characteristic timescale T 2 asso- 
ciated with this noise (over longer timescales the bias 
fluctuates over a larger range- see section 3 for more de- 
tails). If < T^rjj, the noisy bias fluctuates rapidly 
compared to the molecular tunneling dynamics, causing 
incoherent tunneling (this is the 'fast diffusion' limit for 
the spin bath 0). This case is illustrated in Fig. 1, which 
shows the way the two qubit levels are affected in time 
by the fluctuating bias. In the opposite extreme of large 
A Q one has a much smaller noise contribution to 7 ^ of 



70 = N eff /nA T 2 



(12) 



where N e f / is the number of bath spins which are active 
(the exact number depends on the particular system- typ- 
ically N e ff ~ O(N), but the exact fraction N e ff/N can 



vary widely from one system to another) . This noise con- 
tribution 7^ <g; 1 (ie., it only weakly affects coherence). 
The analogue of Fig. 1 for this case would show very 
small fluctuations, which hardly affect the dynamics. 
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FIG. 1: We show the effect of a randomly fluctuating envi- 
ronmental noise bias e(t) (black curve) on a tunneling 2- level 
qubit with tunneling matrix element A D . The 2 levels having 
adiabatic energies ±E{t) , with E 2 (t) = A 2 a + e 2 (t) , are shown 
as red & blue curves. The system can only make transitions 
when near "resonance" (ie., when |e(t)| is ~ A or less, the 
regions shown in green). The resulting dynamics of the qubit 
is incoherent in this case of strong noise. 

(b) Precessional Decoherence: The field about which 
the fc-th bath spin precesses changes each time the qubit 
flips, so that the time evolution of the spin bath states 
becomes entangled with that of the qubit. We can visual- 
ize this process by imagining the precessional motion of a 
bath spin in the qubit field (Fig. 2). Integrating out the 
spin bath then gives decoherence in the qubit dynamics. 
If the "operating frequency" A Q is low, ie., A Q -C E a , 
then this "precessional decoherence" contribution 7^ to 
70 is given by 



_L\2 



7^ = ^E(4/^ 

k 



[if »a4 A ° 



(if 4^<4,A ) (13) 



One gets the second result from the first by a duality, 
switching the roles of uA and uj^r in the derivation of the 
first (cf. ref. 0, App. 2B). 

If instead A D 3> E Q , ie., high operating frequency, then 
also A Q 3> u> \ ojj :. The solution of this weak coupling 
problem is 013|: 



7 £ = (E /A ) 2 /2 



(14) 



and this result is clearly important for the regime of co- 
herent qubit dynamics. 
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(c) Topological Decoherence: When the qubit flips, 
it causes a sudden time-dependent perturbation on the 
bath spins, described by the non-diagonal term ©. This 
induces transitions in the bath spin states, and a corre- 
sponding contribution to the entanglement of the bath 
spins and qubit states. Formally this entangles the topo- 
logical Berry phase of the qubit 01 with that of the bath 
spins, in the same way as for precessional decoherence; 
after averaging over bath states the resulting contribu- 
tion 7^ to 70 is 

k 

where |d?fc| = 7r|co>!|/2f2„is assumed to be small (for gen- 
eral coupling see refs. [71 Il7j|). In general this contribu- 
tion is smaller than the precessional decoherence. 




FIG. 2: We show schematically the motion of a satellite spin, 
in the presence of a qubit flipping between 2 different states 
j) and | J.). When the qubit flips, the qubit field acting on 
the fc-th satellite spin rapidly changes, from 7^ to 7^ (or vice- 
versa). Between flips the spin precesses around the qubit field, 
accumulating an extra "precessional" phase. Averaging over 
this phase gives precessional decoherence. The sudden change 
of qubit field also perturbs the satellite spin phase, giving fur- 
ther decoherence (the "topological decoherence" mechanism 

0). 

We now observe, as has been noted before [3^], that the 
spin bath decoherence rate is always higher at low energy 
(small A ), whereas the oscillator bath decoherence rate 
is higher at high energy (large A D ). Thus there will be a 
'coherence window' at intermediate values of A D , where 
70 is small. 

2c: DECOHERENCE CROSSOVERS IN A 
MAGNETIC QUBIT 

At low T the spin Hamiltonian of many large-spin 
nanomagnetic systems (magnetic molecules, rare earth 



ions, or nanomagnetic particles) reduces from that of 
a tunneling spin S to a simple 2-state form H (f) — 
(A t x + t T z ), with the Pauli spin f acting on the 2 low- 
est spin levels E3 % as m our qubit Hamiltonian . 
The spin gap E g to the next levels is typically ~ 5— 10 K, 
and the 2-state picture is valid at energies <C E g . We as- 
sume henceforth an "easy z-axis" nanomagnet; then the 
'bias' energy e D = gfJ.BS z H*. When e G = 0, the split- 
ting A Q between the 2 "qubit" states |— ), |+) (bonding 
and anti-bonding eigenstates of TL {t)) is produced by 
tunneling between 2 potential wells, with each well hav- 
ing a "small oscillation" energy fl Q ; typically fl Q ~ E g , 
The qubit is thus the result of truncating out the higher 
spin states of the nanomagnetic system, which we should 
schematically in Fig. 3. 

Energy 




FIG. 3: Magnetic anisotropy barrier of a small magnetic sys- 
tem, such as a magnetic molecule. We show the eigenstates 
\S,m) of the longitudinal part TLo{S z ) of the spin Hamilto- 
nian 7i (S). An external longitudinal field H^ xt (or an inter- 
nal field Hi nt ), biases the effective potential by an amount 
£. Adding transverse anisotropy terms to the spin Hamilto- 
nian causes tunneling between the states in the figure. The 
'ground state' tunneling amplitude As between states \S,S) 
and 15, —S) is called A in the text. 

We define the states | f),| J.) (eigenstates of f z ) by 
|±) = I t) ± I {). If the total nanomagnetic spin S is not 
too small, these states correspond roughly to semiclassi- 
cal spin coherent states ^lj , having orientations rv a (here 
o = t, I), which depend on both the internal anisotropy 
field of the nanomagnet, and any transverse external field 
. The splitting A D depends sensitively on . 

The intrinsic decoherence in insulating nanomagnets 
comes from entanglement of the nanomagnetic spin wave 
function with that of the nuclear spins and phonons 0] . 
We first see how to write these couplings in the form 
given in section 2(a). The details for the spin bath are 
a slight generalisation of this form, because the nuclear 
spins are not necessarily spin- 1/2 objects. 

The nuclear spins {Ik} couple to the electronic spins 
{sj} in S (where S = J2j s j) y i a hyperfine couplings 

Hh W = A% l 8fl£ (16) 
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whose form we do not specify here. We then define the 
field-dependent quantity 



UJ 



3 



(17) 



where (s™ )" is the expectation value of Sj when S — > 
Sn a . The energy change of Ifc when S flips from Sn^ 

to Sail is then 2/fewj!, ie., there is a diagonal coupling 

T z w|/fe • I/c between the qubit and Ifc, where Ik is a unit 
vector parallel to the hyperfine field on Ifc. This is just 
the coupling specified in 10} in the last section. 

The external transverse field H^- couples to Ifc with 
Zeeman coupling uik^k 'Ifc, where Wfc-mfc = /ijvH^ and 
rhk is a unit vector along H^. This is the same as the 
coupling given in @. Finally, the interactions between 
the spins in the spin bath add a term 



N N 



n NN ({i k }) = LL v&m 



(18) 



k=l k' = l 



whose effect on the qubit will be handled by assum- 
ing that the spin diffusion caused by this weak inter- 
action adds a "noise" term £(t)r z to the static bias e r z . 
These terms taken together define an effective interaction 
Hamiltonian Hns = H (f) + T^(r,Ifc), where 



V = f z 



+ 52 UJ km k -i k (19) 



To H19[l we also add a 'spin-boson' coupling 



H 



S(j> 



C q ^x x q 



(20) 



to the acoustic phonon coordinate x q . This non-diagonal 
term was already discussed above. 

The energy scale over which the nuclear spin bath 
states operate is just the linewidth E a of the entire mul- 
tiplet of nuclear states coupling to each qubit level [18J . 
It is clear that if the energy bias £ in the problem is less 
than E Q , it will be possible for the system to make tun- 
neling transitions without the aid of the phonons, even if 
£ ^> A Q (see Fig. 4). For this problem, with higher spin 
nuclei, one easily finds that 



k 



4 + 1 , | 

(to 

34 



fc J fc) 2 



(21) 



The acoustic phonon energy scale is the Debye energy 
6b- Now in a nanomagnetic system the ratio E o /ks0D 
can be < 10 , suggesting the a very simple tactic for 
suppressing decoherence. If we tune A D so that kB@D S> 
A c ^> E a , then we will be in the "coherence window" 
mentioned above, where the qubit dynamics is too slow 
to disturb most phonons appreciably, but too fast for the 
nuclear spins to react. 



|S,S> 




|S,-S> 



FIG. 4: The 2 qubit levels connected with the zero-field states 
\S, S) and \S, —S) each couple, via the internal hyperfine in- 
teractions, to a a very large number of nuclear spin levels. The 
result is a multiplet associated with each qubit level, which 
usually has a Gaussian density of levels, with a half-width E 
(see text). 



To substantiate this idea, we generalise the low field 
(A < E ) calculations of nanomagnetic dynamics [l8j |. 
where incoherent tunneling relaxation is found, to the 
high-field regime A Q E a . Because the {^f.} <C the 
Zeeman couplings {w^}, and <C A Q , this dynamics 
can be solved 0, EH, by expanding O in wjj/A. If 
we first ignore the noise term, we get a new effective 
Hamiltonian 



Hnr — 



A + 



kk' 



2A 



— (Ik ' Ifc)(4' • Ifc' 



+ V^mfc.Ifc + 0((^) 4 /A^) 



fc 



(22) 



We can easily generalise the derivation of the result given 
in (|14fl for 7^ to the case where the nuclear spins have a 
spin modulus Ik > 1/2. One finds easily that 




(23) 



Thus we see that the result <|14[) for precessional deco- 
herence is generally valid, regardless of the size of the 
nuclear spins. 

We now look at the 2 other contributions to the de- 
coherence rate. The nuclear spin transitions induced 
directly by electronic spin flips add a contribution 7^ 
to 70. However when A ^> E , the ratio 7^/7^ ~ 
0(&? /£l 2 ) <C 1, ie., the precessional decoherence always 
dominates over the topological decoherence. 

The third contribution 7^ comes from nuclear spin 
noise (the term £(t)r z ). When A D ^> E Q , these fluc- 
tuations are extremely slow compared to A Q ; typically 
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Ti ~ msecs at low T, where T2 is the typical transverse 
nuclear relaxation time for the N nuclei controlling these 
fluctuations. One then gets 7^ = N/nA T 2 ; we will see 
this is very small. 

This summarizes the nuclear spin terms. The phonon 
contribution is as described above in ©, and in the low 
T limit we are interested in, where k B T < A Q , we have: 




[(^ A ) 2 /e^] (24) 



At temperatures above A D the phonon decoherence rate 
increases. 

Now, since 7^ dominates nuclear spin decoherence, we 
can get an estimate for the optimal decoherence rate 7™" 

by simply minimizing 7^ + 7^ with respect to A D , as- 
suming k B T < A , to get: 

7^ m « V2Sn E /6 2 D (25) 

at an optimal tunneling splitting A° pt : 

Ai°rt „ eniEo/V^SSlo) 1 / 2 . (26) 

We see that decoherence is optimised for a given S by 
making E and fl small, and Ojj large, within the con- 
straint that fi > A c > k B T. If k B T > A Q we get a 
different (less favorable) answer. 

The detailed application of this kind of result to a 
nanomagnet, in cases where one knows something about 
the couplings, is in principle very useful for designing 
magnetic qubits. The tunneling splitting is most easily 
modified just by applying a magnetic field transverse to 
the easy axis- this can be used to tune A D over many or- 
ders of magnitude (see Fig. 5) . For more details of such 
applications, see ref. [25| . 

1E+1 -_ 




FIG. 5: The tunneling splitting |A | in the Fe-8 molecule, 
which is a good example of a tunneling nanomagnet. The 
tunneling anisotropy potential for the spin is biaxial, with 
an easy z-axis and a hard :r-axis. The tunneling splitting is 
shown as a function of a transverse field H± oriented in the 
a;y-plane, at an angle <f) from the hard axis. 



3: INTERLUDE-SPIN BATH DYNAMICS 

In the results given above for decoherence rates, the 
intrinsic dynamics of the oscillator and spin baths played 
only a secondary role. In both cases it was assumed that 
phase information exchanged between system and bath 
was lost once it was taken up by the bath. This may not 
always be realistic, particularly in the spin bath- one can 
easily imagine situations in which the spin bath is cycled 
so as to recover some of this phase information, and in 
NMR this is actually done (eg., in 'multiple quantum 
coherence' experiments 42] ) . Even for oscillator baths it 
is well known that anharmonic oscillator couplings can 
allow the bath to hold information in certain modes for 
very long times- now a very well-studied phenomenon 

The question of the intrinsic bath dynamics is also im- 
portant when one looks at the crossover between quan- 
tum and classical relaxation (next section). Therefore 
here we clarify what is and is not contained in the mod- 
els we use. 

As described by equation J2J, the oscillator bath has 
the very simple dynamics of N independent oscillators, 
with frequencies {tu q }. Coupling to the 'central sys- 
tem' hardly changes this dynamics, since the coupling 
strengths c q ~ 0{N~ 1 / 2 ). Thus nowhere in this model 
does energy and phase relaxation occurs in the bath- 
both simply accumulate independently in each mode. 
In reality anharmonic couplings cause rapid relaxation 
in systems of extended modes like electrons, phonons, 
magnons, etc; the only system to which © strictly ap- 
plies is a bath of photons in a vacuum, in which the very 
weak vacuum polarisation-induced photon-photon inter- 
actions have been dropped. The reason that models like 
the spin-boson model work in most (but not all!) cases is 
just because this relaxation is usually fast- once energy 
or phase information has gone from the central system 
into a particular mode, it is rapidly diffused into other 
modes and so hard to recover. 

In the case of the spin bath dynamics one has to be 
more careful. As discussed in detail in ref. [7J, one 
can classify the large number of spin bath states into 
'polarisation groups', defined by their total polarisation 
M = Zk ■ &k along some set of axes where the 

{iifc} are unit vectors. If the 'external' field strengths 
{co^} are weak, so that u^r -C then it makes sense 
to have these axes along the direction of the local fields 
acting on the {crfc}, coming from the qubit, ie., to make 
Zk = Ik (compare eqtn (0}). If on the other hand the ex- 
ternal fields dominate, one instead assumes that Zk = rhk, 
ie., the axes of quantization defining M are just the ex- 
ternal field directions. If these fields hj, = uj^rhk are 
indeed defined by some strong external magnetic field 
H = h H , then Zk — h a for all spins in the bath, ie., 
the axis defining the bath polarisation is just the exter- 
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nal field direction. However we emphasize that in many 
cases the fields hfc may have nothing to do with any ex- 
ternal magnetic field. For example, the {cr^} might refer 
to a set of defects in a glass- in this case they couple to a 
strain field like electric dipoles. Even in the case of real 
spins the hfc may not be external fields but internal ones. 
For example, many tunneling nanomagnets do not tun- 
nel between oppositely spin states, because the magnetic 
anisotropy field does not have uniaxial symmetry. In this 
case the hyperfine field on the nuclear spins does not flip 
through 180° when the nanomagnet tunnels. One then 
resolves the hyperfine field into 2 components; hfc defines 
that component which does not change during the flip, 
while ujUk defines a component which flips through 180°. 

The point of defining polarisation groups in this way 
is that it then makes sense to define 2 relaxation times 
Ti and T 2 , such that T\ defines the relaxation of M, and 
T2 the spin bath relaxation within a given polarisation 
group. Note that even if the spin bath is actually made 
up of nuclear spins, these times are not the Ti,T 2 times 
measured in a typical NMR experiment, which only looks 
at a single nuclear species at a time- the times here refer 
to the whole spin bath. One can imagine an NMR exper- 
iment which polarizes a particular nuclear species, and 
then observes rapid relaxation of this polarisation into 
other nuclear species, via the inter-spain interactions- 
even though the total spin along the field is conserved 
(so that the spin bath T\ defined here is still infinite). 

Usually one expects T\ ^> T 2 at low T, because the 
interspin couplings V^a^a 1 ^, can mediate transverse 
spin-spin relaxation processes, ie., cause spin diffusion 
in the spin bath and contribute to T 2 , whereas T\ pro- 
cesses typically require interaction with some external 
system. The most typical case where the bath spins 
inter-communicate by magnetic or electric dipolar inter- 
actions is actually hard to analyse theoretically, because 
these interactions are marginal in 3 dimensions (integrals 
of the form n?~ J ^ 3r j/ r j appear in the calculation of 
the relaxation dynamics of n coupled bath spins, so that 
multi-spin couplings are just as important as pairwise in- 
teractions, and distant spins as important as nearby ones 
Thus we would not usually try to calculate the spin 
bath f 2 . 

Let us now emphasize one of the crucial differences be- 
tween the spin and oscillator baths. This is that because 
the {uA } are not weak, the spin bath dynamics depends 
very strongly on what the central system is doing. In 
the case of a central qubit one can imagine 2 extreme 
scenarios: 

(i) We freeze the qubit dynamics by applying a longi- 
tudinal field e c 3> A D . Then the spin bath dynamics is 
described by the effective Hamiltonian 

Use = £ b fc • <r fc + £ V$o%4 (27) 

k kk' 



where the static fields {bfc} are given by 

b fe = LUjtm k ± upk (28) 

with the sign ± depending on whether the central qubit 
is frozen in the | T) or | J.) state. 

We see that in this extreme case the bath would have 
its total polarisation conserved, provided we defined the 
polarisation groups using axes {zk} parallel to the {bfc}. 
Actually we would not normally do this, but one can 
easily imagine a situation in which the applied field is 
either very strong or very weak, and then M would be 
almost exactly conserved, ie., T\ would be very long. 

Now suppose we switch on the qubit dynamics- the 
easiest way to do this is to remove the bias field e . Then 
the nuclear bath finds itself subject to a quite different 
time-dependent Hamiltonian, of the form 

Ksb® = E^mfc-afc + Y. V ^^i4 

k kk' 

+ T z (t)J2"Pk-Vk (29) 

k 

where the timc-dcpcndcnt variable r z (t) is jumping back 
and forth between ±1 with some correlation time ~ 
1/A D . This Hamiltonian described the spin bath now 
subject to an external 'telegraph noise', which causes 
transitions between different polarisation groups of the 
system- this will be true no matter how these polarisation 
groups are defined, provided both the {uj^:} and the {wf } 
are non-zero. Suppose, for example, that we have a weak 
external field. Then M is defined as M = J^k ™ k ' as 
discussed above. We can imagine an initial state where 
all bath spins are oriented parallel or antiparallel to the 
initial local fields- but as soon as the qubit flips, they 
begin to precess. As discussed above, this is what causes 
precessional decoherence. 

However, as observed already above, the motion of the 
central qubit depends itself on the spin bath dynamics- 
we must never forget that the telegraph noise acting on 
the spin bath depends in turn on the spin bath state. The 
simplifying feature is that energy conservation imposes 
a simple constraint on the allowed bath dynamics. In 
general the qubit will be off resonance by some energy £, 
the sum of an external field contribution and the internal 
field from the spin bath: 

t = e + J2"l l k-Vk (30) 
k 

In the simplest case where the couplings u>) are either 

dominated by a single value uio , or else all cluster around 
this value (this happens in , eg ., rare earth magnets like 
the LiHo x Yi^ x Fi system |52j . where the Ho hyperfine 
coupling to the Ho nuclear spin is much larger than its 
coupling to the other nuclear spins), we have approxi- 
mately that £ ~ e Q + lo\m (with some spread around the 
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value Wo, caused by dispersion in the values of the w\). 
Now for the qubit to make transitions we require that 
the initial and final energies ±£ (for r z = ±) be the same 
within A . This means that if e ^> A , the central qubit 
can only flip if the spin bath absorbs the extra energy. 
However this will in general involve a change of 2M in 
the net bath polarisation, such that cj|J|M| ~ |e G |, each 
time the system flips. Thus at least M bath spins have to 
flip- the time-varying field of the qubit must drive these 
transitions. 

A formal calculation of the spin bath dynamics incor- 
porates this constraint using a projection operator 0, 
involving a dummy variable £: 



n M = 



JV 

fc=i 



' Jo 2TT 



-M) 



(31) 



which restricts all bath states to the M-th polarisation 
group. Suppose now that the spin bath starts off with 
polarisation M = M a , and that to maintain resonance, 
the net polarisation must change by ±2M each time the 
qubit flips from | |) to | J.) or vice- versa, ie., it cycles 
between M a < — > M a —2M. The simplest case arises if we 
ignore the interaction between bath spins completely, ie., 
let — > in 129( 1 . Then one can write the dynamics 
of the bath in terms of operators T n and acting on 
the spin bath in the presence of n flips of the qubit. The 
T n are given by 



T„ 



>U 



u 



(32) 



(involving a set of n dummy variables) and the {Uk} de- 
fine the change in the wave-function of the bath spins 
caused by the sudden flip of the qubit from one orien- 
tation to another, ie., the mismatch between in and out 
states: 



wr t }) = f[Uk\{vi n }) 

k=l 



(33) 



Suppose, eg., that ujjt <C for all bath spins, ie., the 
field on each bath spin almost exactly reverses during 
each flip, through an angle 180° — 2(3k- Then Uk is just 



(34) 



Suppose we now want to write down the amplitude 
for the spin bath to start in the polarisation group M 
and finish in the same polarisation group. This can only 
happen if the qubit flips 2n times. The amplitude is then 
the sum of a term 

\2n 2n 



-q I d& c _iM ($ 2 „+$ 2 „- 1 +...+gi) 



x e 



(iA ty 

(2n)! 



-J- 2 



(35) 



acting on the initial state of the spin bath, in which the 
qubit is assume to start and finish in the same state | f), 
and another terms in which it starts and finishes in the 
state | i). 

The result of form (|35|l is only complete if energy con- 
servation requires that the polarisation change by ±M 
each time. However in most cases there will be a wide 
range of values of w|, rather than a single dominant value, 
and so each polarisation group will be widely spread in 
energy space, and a large number of polarisation groups 
will have states at a given energy £ (ie., the groups will 
strongly overlap in energy space). In this more general 
case we should sum over transition amongst these groups 
with the appropriate weighting- the details are an obvi- 
ous extension of what has just been described. In this 
way we can give a theoretical evaluation of the time T\ . 

To calculate the full dynamics of the spin bath we must 
also include the action of the interspin interaction V k a y ■ 
This enables transitions amongst the different bath states 
inside the same polarisation group, even when M does 
not change, ie., an evaluation of the spin bath T^. We do 
not discuss here how such calculations may be done. 

We underline here again the most important point of 
this interlude- that the spin bath controls the qubit dy- 
namics, deciding whether the qubit may flip or not- but 
in its turn the qubit drives the spin bath dynamics. 



4: QUANTUM TO CLASSICAL CROSSOVER 

As noted in the introduction, a great deal is known 
about how the dynamics of a single quantum system, 
coupled to a thermal bath, changes as one raises the bath 
temperature 0,0. Many analyses use an oscillator model 
to describe the bath. However the models usually used 
to discuss this problem are restricted in certain ways. It 
is assumed that the bath stays in equilibrium during the 
times of interest, so that the internal relaxation times 
in the bath must be short compared to the timescale 
relevant to the dynamics of the quantum system, and 
energy given by the quantum system to the bath rapidly 
moves away, redistributing itself amongst bath modes. 

The problem with such models is that at low T the ba- 
sic assumption of short internal relaxation times breaks 
down. This is of course well known and has been studied 
theoretically in, eg., spin glasses 0,0], dipolar glasses 
[lif , some models of low-T phase nucleation > an d in 
the various relaxation bottlenecks existing in magnetic 
systems 0] (to name only a few examples). However 
the problem is really generic to low-T physics (a feature 
well-known to low-T experimentalists |49j, because it is 
the main obstacle to cooling to very low T). In fact, in al- 
most all systems apart from pure liquid 3 He and i He, the 
thermal bath has most of its low T energy and entropy 
locked up in localised excitations. The relaxation of en- 
ergy and entropy in and out of these modes ranges from 
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/zsecs to centuries. From a theoretical standpoint these 
facts are not surprising- they often arise in systems of 
local modes when couplings and fields are random (par- 
ticularly when there is frustration of some kind in the 
interactions). 

What this means here is that we should reconsider the 
whole problem of the quantum-classical crossover, using 
models which have these localised modes built in to them 
from the start. Consider, eg., the standard 2-well system, 
coupled now to both oscillator and spin baths. The os- 
cillators represent delocalised modes, which at low T are 
few in number but can move energy around quickly, and 
the spin bath represents the localised modes, which con- 
tain all the energy and entropy. A toy model for such a 
system is 

H = H {P : Q)+H° SC + H S B + V int 
H (P,Q) = P 2 /2M + U (Q) 

Vint = Q^2,c q x q + ^2 wf4 • <7fc] (36) 

q k 

in which the oscillator and spin baths are as before (see 
eqtn. ©), and the couplings are simple diagonal ones to 
the coordinate of the particle. One can have more com- 
plicated coup ling s, and in general one should also add 
counterterms |5fl | to Ti to renormalise the 2-well tunnel- 
ing potential back to U(Q). 

Already in this toy model one begins to see how the 
quantum-classical crossover will work. At low T the dy- 
namics of the particle will be governed by its coupling to 
the spin bath, in the way previously described. Raising 
the bath temperature has little effect on the spin bath dy- 
namics unless the {wj!} are very large- they are already at 
high T compared to their basic energy scales {wj(,a->^}. 
However eventually the effect of thermal transitions of 
the oscillators to higher levels of the 2-well system be- 
gins to take effect. In the absence of the spin bath one 
sees a crossover to activated behaviour around a temper- 
ature T Q ~ cj /27t, where lo q is the (renormalised) small 
oscillation frequency of the system in one or other of the 
wells PJ . The width AT of the crossover depends on the 
details of the potential, the bath coupling, etc., but it 
will not be less than ~ 0(T /n), where n is the number 
of levels below the barrier. 

However the spin bath introduces new timescales in 
the problem, viz., the T% and Ti introduced in the last 
section. Now these timescales, depending as they do on 
the dynamics of the central system itself, will decrease 
rapidly as we raise the temperature- the more rapid fluc- 
tuations of the central system, caused by coupling to 
thermally excited oscillators, stimulate more rapid tran- 
sitions in the spin bath. This indirect effect of the os- 
cillator bath on the spin bath dynamics, acting through 
the central system, is of course well known in NMR. In 
any case, we see the spin bath can compete with the os- 
cillator bath over a rather wide range of temperatures in 



controlling the dynamics of the central system. This tells 
us that we may expect a much wider crossover between 
quantum and classical behaviour than occurs when one 
only deals with an oscillator bath environment (or only 
a spin bath environment). 

Rather than give a general study of this crossover here, 
which is rather lengthy, we now present instead some rel- 
evant results for magnetic molecular spin relaxation, for 
which there also exist fairly detailed experimental results. 
These results actually capture some of the more general 
features of the problem. 

4A: QUANTUM-CLASSICAL CROSSOVER FOR 
MAGNETIC MOLECULES 

It has been understood for many years that the ther- 
mally activated spin dynamics of insulating magnetic 
ions and nanomagnets (including large spin magnetic 
molecules) is driven by coupling to the phonon bath. Ex- 
perimental investigations in the last 10 years of the tun- 
neling relaxation dynamics of various magnetic molecules 
has led to many theoretical attempts to understand 
the temperature dependence of this behaviour in terms 
of the spin-phonon coupling to single tunneling spins 
[53l IHil l55l I56L |H3 • Such calculations give a typical 
crossover from some straightforward low-T tunneling to 
higher-T activation- the details tend to be be messy be- 
cause of the presence of many levels and different spin- 
phonon coupling. These calculations must clearly ap- 
ply at some sufficiently high temperature that neither 
nuclear spins nor intermolecular dipolar interactions are 
relevant. However they are not directly applicable to the 
experiments in the crossover regime, nor in the quantum 
regime. This is clear both on theoretical and experimen- 
tal grounds, as follows: 

From a theoretical point of view the interplay be- 
tween nuclear and phonon couplings on a single nano- 
magnet cannot be ignored, especially given that the nu- 
clear dynamics is fast. In experiments the intermolec- 
ular dipole interactions couple the relaxation of differ- 
ent molecules, so that they can only relax independently 
when ksT ^> Vd, where Vd is the strength of these dipo- 
lar interactions. Thus until we reach this rather high 
temperature (which in most experiments is well above 
the temperature T a defined earlier), both intermolecular 
dipole and hyperfine coupling to the nuclear spins must 
be included on an equal footing with the spin-phonon 
interactions. 

From the experimental point of view the need for this 
is obvious. Even well above T a the relaxation is non- 
exponential in time (as in the quantum regime, but now 
with T-dependent characteristics) , showing the molecules 
do not relax independently. Moreover, the 'hole-digging' 
phenomenon in the distribution of internal fields, caused 
by nuclear spins, also survives at temperatures well above 
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T Q (although with a time-development that becomes 
rapidly T-dependent), demonstrating that the nuclear 
spins are still partially controlling the tunneling dynam- 
ics. 

To begin analysis of this problem we note that be- 
cause the nuclear spin bath and the phonon bath do 
not interact with each other directly (but only via the 
central molecular spin), we can treat their relaxation 
rates as independent. We assume that each nanomag- 
net has a spin Hamiltonian of easy axis form H a (S) — 
H*(S z ) + Vj-(S±), so that the transverse term Vj~ causes 
tunneling between the eigenstates \m) of the longitudinal 
part (defined by H%\m) — Em\m)). Let us also assume 
that the applied longitudinal field is small, so that levels 
\m) and | — m) are near resonance (and level |m) is not 
near resonance with any other levels). Now suppose we 
start with an ensemble of such nanomagnets, all of them 
having m = S, ie., so the system is completely polarised. 
If we ignore any coherence effects (ie., assume incoherent 
relaxation), and also ignore intermolecular interactions, 
then the kinetic equation for the system reduces to 

P£>(£,r,t) = - £[r m , m (£,T)P«(£,r,t) 



r w (£,T)P«(£,r,f)] (37) 



(r_ m , m = r TOi _ TO = T m ) where Pfh'(£,r,t) is the 1- 
molecule probability distribution, describing the proba- 
bility to find a molecule at position r, in state |to), in a 
bias field £, at time t; and r mm /(£, T) is the rate at which 
nanomagnets in a local field £ make transitions from state 
|m') to |m), under the influence of both phonons (at a 
temperature T) and nuclear spins. The assumption of 
non-interacting phonon and nuclear baths implies we can 
write: 

Iw (£, T) = r£ m , (i, T) + Ti m , (£, T) (38) 

with individually defined nuclear spin- and phonon- 
mediated relaxation rates. The system can move up or 
down levels on the same side of the barrier by emission 
of phonons. Of particular interest here are the inelastic 
tunneling rates out of level |m) in a bias field, to the 
other side of the barrier. The phonon-mediated process 
of this kind has the form 1531: 



r£(0 



A? n W m (T) 



(39) 



where W m (T) is that part of the linewidth of the m- 
th level caused by phonon-mediated intra-well processes, 
£m = g/ismiP is the bias energy between levels | ± m) 
(with the field H z the sum of internal and external fields), 
and A m is the tunneling matrix element for tunneling 
between these same levels. 



On the other hand for the nuclear spin bath-mediated 
rate we will use here the temperature- independent form: 



C(0 



(m) 



(m) 



(40) 



where E^" 1 ^ is a generalisation of the quantity E a which 
plays a role in nuclear spin-mediated tunneling in the 
quantum regime- roughly one has Eo ~ (m/S)E (cf., 
cqtn. pi|l). This quantity measures the range over which 
the nuclear spin bias is being swept, either by internal 
spin diffusion or by transitions caused by the nanomag- 
netic dynamics (the mechanism described in the last sec- 
tion). The factor G^ l) = exp{-AlJ2(E { m) ) 2 } follows 
from the Gaussian spread of the nuclear multiplet, and 
simply says that (£ ) vanishes when A m becomes large 
in comparison with Eo^ ■ 

The rationale for using (|40|1 is that we are interested 
here in time scales long compared to the T\ and T2 dis- 
cussed in the last section. In this case we expect that 
the system is able to cover the whole range of states in 
the nuclear spin manifold surrounding each level, and so 
we can simply weight these according to their number, 
ie., according to a density of states. Clearly this approx- 
imation breaks down if we are interested in shorter time 
scales- we will not consider this problem here. We can 
in fact go further- since at time scales longer than all re- 
laxation times, the nuclear bias will fluctuate across the 
whole range of nuclear states, of width -Eg™-*, we can in- 
corporate this into the phonon rate as an average- for 
example, when E^ n ^ >> max{A m , KW m }, the phonon 
rate becomes 



E 



(m) 



77 -e/2(^ m) ) 2 



(41) 

Actually this result turns out to be valid even when the 
spin bath dynamics is slow compared to the timescale we 
are interested in- but with one simple modification. The 
point is that even if a given molecule in a dipolar bias 
field £d, coming from the other molecules, cannot find 
resonance, because the spin bath brings it too slowly to 
resonance, nevertheless some fraction x of the molecules 
in a field £ will find themselves in a compensating nuclear 
bias field near — £d, near enough so that the molecule is 
quickly brought to resonance. 

We see that in this very simple approximation the spin 
and oscillator bath-mediated relaxation processes already 
influence each other strongly, albeit in a rather trivial 
way. Because l|39|l reduces the contribution of the higher 
levels to the relaxation, which in the usual theory of the 
crossover take over very quickly from the lowest levels as 
one goes through T = T , the net effect is to broaden the 
crossover. 

To illustrate this it is useful to show results for a par- 
ticular system. We again choose the Fe-8 molecule, al- 
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ready discussed above in the context of low-T decoher- 
ence. In Fig. 6 we show the contributions from the differ- 
ent levels to the relaxation as a function of T for the Fe- 
8 system. To give results comparable with experiment 
we have generalised the kinetic equation l|57|) to include 
dipolar interactions- this development is a rather messy 
but fairly straightforward adaptation of the method used 
m ref. 591. The main effect of adding these interactions 
is however not to change the width of the crossover, but 
rather to change the time dependence of the relaxation- 
this it is fairly complex, even in the quantum regime |59| , 
and not relevant to the present study. 



x3< m >x(1-M(t))/M(t) 



6 mK, M(0)=1 



5(x=0.1) 
6(x=0.1) 

7(Hl=0) 




FIG. 6: The relaxation rates as a function of temperature 
for relevant values of S z — m in the case of the Fes system. 
The two curves for m = 7 correspond to different transverse 
fields, viz., (i) H± = (filled triangles); and (ii) Hj_ — 55 mT 
(open flipped triangles) . The inset enlarges the m = 10 curve 
around T >T , where T is the temperature at which phonon 
transitions would normally cause the tunneling dynamics to 
rapidly crossover to activated behaviour. 

The most important point to be noted from Fig. 6 
is that the crossover is now very wide- it extends from 
T ~ 0.4 K up to T ~ 1.7 K, above which temperature 
the contributions from level m = 5 begin to dominate. 
It will be surprising to those familiar with the standard 
theory of the quantum-classical crossover that it is an 
intermediate level that dominates. The crossover is so 
wide for 2 reasons, viz., 

(i) The reason noted above, ie., the spreading of the 
levels by the nuclear spin hyperhne coupling (and in the 
case of this calculation, also the dipolar fields, so that 
now the spread is even greater than E 1 ^). This em- 
phasizes the role of the lower levels, more than would 
otherwise be the case; and 

(ii) We see already that in the basic phonon transition 
rate in (|39J) there is a saturation in the rate for the high 
levels (having small m and large A m ). This is because 
the typical bias £ m in this formula, and in the more gen- 
eral formula l|41|) , will now be either a dipolar or nuclear 



hyperhne bias, which is much larger than the phonon 
linewidth W m . Once A m exceeds this bias, the presence 
of A m in the denominator stops the rapid increase of the 
rate- this happens before one reaches the very highest lev- 
els, and is the basic reason why intermediate levels domi- 
nate the relaxation over a wide temperature range. In the 
case of Fe-8, this happens for m = 5, but clearly could 
happen for some intermediate level in a different system. 
In most of the earlier papers on phonon-mediated tun- 
neling relaxation of nanomagnets, the factor of in 
the denominator was not included- this led to a quite 
different picture of the crossover. 

Actually experiments in these systems do show a very 
wide crossover. The detailed comparison between theory 
and experiment is rather interesting, since one may anal- 
yse both the relaxation as a function of time (and how 
the form changes with T, along with rate) and also the 
T-dependence of the hole-digging dynamics j^jj . 



SUMMARY 

This paper has not attempted a complete study- in- 
stead we have tried to make some general points about 
2 kinds of crossover, uncluttered by too much detail for 
particular systems. To illustrate the general remarks it 
has been nevertheless useful to give results for the Fe-8 
molecule. It goes without saying that the detailed calcu- 
lations for this system (and others, such as SQUIDs) are 
rather lengthy, and including them would have obscured 
the points we wishes to make. 

The 2 main points are that 

(i) The presence of the spin bath is very bad for de- 
coherence when the basic energy scale A Q of a qubit is 
small. This means that the spin-boson model is not ap- 
plicable at all in this regime- one must use a 'central spin' 
model 0- However if we increase A Q , the spin bath deco- 
herence effects fall off extremely rapidly, and eventually 
become negligible. At this point the spin-boson model 
becomes applicable, and decoherence begins to increase 
again as one further increases A D , a well-known feature 
of the model Q] . 

(ii) In the standard quantum-classical crossover that 
occurs as one raises the bath temperature for a tunnel- 
ing system, the presence of the spin bath is again very 
important. The main effect it has is first to completely 
change behaviour in the quantum regime, and then to 
enormously broaden the usual rather sharp crossover that 
exists when one only deals with an oscillator bath. 

It is quite clear that these results are only the begin- 
ning of a proper study of the way in which spin and 
oscillator baths work together. The remarkable edifice of 
theoretical work that has been constructed around the 
spin-boson and related models should be a very nice 
model for what interesting paths remain to be explored 
in this area. 
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